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Abstract 

Damping rates of multipolar, localized surface plasmons (SP) of gold and silver nanospheres 

of radii up to lOOOnm were found with the tools of classical electrodynamics. The significant 

increase in damping rates followed by noteworthy decrease for larger particles takes place along 

with substantial red-shift of plasmon resonance frequencies as a function of particle size. We 

also introduced interface damping into our modeling, which substantially modifies the plasmon 
P-i 

damping rates of smaller particles. We demonstrate unexpected reduction of the multipolar SP 

damping rates in certain size ranges. This effect can be explained by the suppression of the 

nonradiative decay channel as a result of the lost competition with the radiative channel. We 

\ show that experimental dipole damping rates [H. Baida, et al., Nano Lett. 9(10) (2009) 3463, 
o t 

and C. Sonnichsen, et al., Phys. Rev. Lett. 88 (2002) 077402], and the resulting resonance 
quality factors can be described in a consistent and straightforward way within our modeling 
extended to particle sizes still unavailable experimentally. 

Keywords: plasmon damping, radiation damping, interface damping, surface plasmon res- 
onance, multipolar plasmons, multipolar plasmon modes, Mie theory, gold nanoparticles, sil- 
ver nanoparticles, nanosphere, nanoantenna, nanophotonics, plasmonics, size dependent optical 
properties, SERS, SP. 

1. Introduction 

The response of noble-metal nanoparticles to electromagnetic (EM) excitation is dominated 
by resonant excitations of localized surface plasmons (SPs) (see: [l|, 0, [si, Q, 0] for reviews). When 
a metal particle is illuminated at corresponding resonance frequency, notably strong surface- 
confined optical fields can be generated. This property is applied in surface-enhanced Raman 
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high-resolution microscopy [9] or improvement of plasmonic solar 
cells |lOl . I 111 . Il2l [13L [15I ] . Excitation of SPs at optical f req uencies , non- diffract ion-limited 
guiding of them (e.g. via a linear chain of gold nanospheres 

0,00, 0, 0, 0] ) and transferring 
them back into freely propagating light are processes of great importance in many applications. 
Gold and silver nanoparticles are important components for subwavelength integrated optics 
and high- sensitivity biosensors 0, 0, 0] due to their chemical inertness and unique optical 
properties in the visible to near-infrared spectral range. 

Metal nanospheres are the simplest and the most fundamental structures for studying the 
basis of plasmon phenomena. Understanding the resonant interaction of light with plasmonic 
nanoparticles is also essential for applications, e.g. for designing useful photonic devices. The 
modeling of plasmon properties under the simplest electrostatic (quasistatic) approximation is 
valid only for particles much smaller than the wavelength of light; than size dependence of 
plasmon properties can be neglected. If a nanosphere is larger, the properties of the surface 



plasmons supported by this structure become dependent on the size of the particle 
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27|. Size dependence of the retarded electronic response of a plasmonic particle is 



determined by the properties of metal at the SP resonance wavelength and the presence of the 
metal-dielectric interface. 

Nanoparticles much smaller than optical wavelength first of all exhibit dipolar surface plas- 
mon oscillations. Higher order multipolar resonances appear as new features in the optical 
spectra, at frequencies higher than that of the dipolar plasmon frequency. The spectral sig- 
natures of these higher order plasmon resonances were observed in elongated gold and silver 
nanoparticles 
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30] first. Fewer experimental investigations directly demonstrated the 
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multipolar character of the observed resonances in spherical particles 

Controlling the spectral properties of a plasmonic nanosphere for technological or diagnostic 
applications is not possible without knowing the direct dependence of plasmon properties on 
particle size. Still it is believed (e.g. js|,0]) that existing theories do not allow the rigorous direct 
calculation of the frequencies of the multipolar SP modes. However, apart from the plasmon 
resonance frequencies, also the plasmon decay rates (damping times) are essential in controlling 
the spectral response of the plasmonic particle. Total damping rates of multipolar SPs |27| define 
absorbing and emitting properties of a plasmonic nanoantennas which can be tuned by particle 
size. Nanoantenna of the right size can serve as an effective transitional structure which is able 
to absorb or to transmit electromagnetic radiation in plasmonic mechanism. The knowledge 
and modeling of plasmon damping effect as a function of particle size are thus of key interest for 
the development of plasmonic nanosystems. The central feature here is the actual enhancement 
of the electromagnetic field. It served as a stimulus for us for developing a strict and direct 
size characterization of multipolar plasmons resonance frequencies and plasmon damping rates. 
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Our extended modeling up to the radius R = lOOOnm, and up to SP multipolarity I = 10 
allows us to determine the intrinsic properties of plasmonic spheres in the size range which has 
never been studied before. We do not apply any restrictions on the particle radius in relation 
to the wavelength of light; our study goes far beyond the particle size for which the quasistatic 
approximation is justified (e.g. 0, 0]). 

The main reason of size sensitivity of SP which we discuss (usually referred to as the retar- 



dation effect ([33j and references therein)) is responsible for important red-shift of the plasmon 
resonance frequency for larger nanoparticles. It is accompanied by the significant increase in 
damping rate followed by noteworthy decrease for larger particles. In addition, we discuss the 
effect called the surface electron scattering effect [lj] or interface damping 34 1 which causes the 
substantial modification of the plasmon damping rates in particles of sizes comparable or smaller 
than the mean electron free path. Sometimes called the "intrinsic size effect", it is describe d by 
the 1/R dependence added to the electron relaxation rate 7 of the dielectric function [l|, 5], |35| . 
We perform both types of modeling including or neglecting the effect of surface scattering, in 
order to examine the causes of the size-dependent plasmon features. 

Our extended modeling for large particles reveals new features of the total plasmon damping 
rate. We demonstrated, for the first time in spherical particles, the effect of reducing of mul- 
tipolar SP damping rates below its low size limit which is equal to the nonradiative damping 
rate. Suppression of the total damping rates proven to exist in certain size ranges allowed a new 
insight into the role of radiation damping in the plasmon decay mechanism. We also describe 
changes of the quality Q-factor of SP multipolar resonances with particle size. Q-factor is used 
for determination of the local field enhancement [34J and effective susceptibility in nonlinear 
optical processes j^ . 

We also confront the size characteristics resulting 



results obtained for the dipole plasmon for gold |32|, |34j and silver [37| particles of different sizes 



iom our modeling with the experimental 



We show that the measured dipole damping rates and quality factors |32|, [34 1 can be described 
in a consistent and straightforward way within our modeling which deliver also predictions for 
particle of size range still unavailable experimentally. 



2. Size dependence of SP parameters derived from optical spectra 

It is widely accepted, that excitation of surface plasmons (collective free-electron oscilla- 
tions) is responsible for the commonly observed pronounced peaks in the optical absorption or 
extinction spectra. The spectra collected for nanoparticles of various sizes are used as a source 
of data allowing to reconstruct the change of the resonance (usually dipole only) position with 

The width of the peak is related to the (dipole) plasmon damping rate 



size (e.g. 
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Mie theory delivers an indispensable formalism enabling the description of scattering of a 
plane monochromatic wave by a homogeneous sphere of known radius, surround by a homoge- 
neous medium 
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45j . Mie formulas allow to predict the intensities of light scattered in 
a given direction or the absorption and scattering cross-section spectra for particles of a chosen 
size. The Mie spectra can be compared with the spectra (intensities of absorbed or scattered 
light) measured experimentally for particles of the same size. 
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Figure 1: Absorption C a f, s (u;) and scattering C sca t(uj) cross-sections for nanospheres of radii a) 8nm and b) 80nm 
(Mie theory). Figures illustrate suppression of absorption in larger particles. Broadening, deformation and shift 
of the maxima in the scattering spectra of larger particles in comparison with the SP resonance eigenfrequencies 



ui[ is also shown. 



However, predicting the size dependence of multipolar plasmon resonance frequencies, Mie 
scattering theory is an inconvenient tool. SP size dependence can be determined only indirectly, 
by laborious derivation of positions of maxima in the consecutive spectra collected for various 
particle radii. The same applies to derivation of the size dependence of the (dipole) SP frequency 
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corresponding to the peak position from experimental data. In addition, the maximum in 
the spectrum ascribed to SP resonance can be broadened, deformed and shifted in respect 



to the frequency of plasmon oscillation mode 26|, |27|, as illustrated in Figure [TJ That can 



affect the plasmon resonance position obtained from the experimental or calculated scattering 
spectra. It is even more difficult to determine the SP damping times from the width of the SP 



resonances, manifested in th e optica 
size dependence then (e.g. 
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spectra of larger particles, and to understand the derived 
3810. The processes leading to plasmon damping are the 



subject of extensive research and debate (e.g. [l|, 34, 0, 40, 0, 47, 48] . The homogeneous 
linewidth of the SP resonance was connected to the SP damping time. The properties of the 
SP are strongly influenced by this parameter. Mie scattering theory is not sufficient to define 
a general rule describing the size dependence of SP parameters; the damping rates and the 
resonance frequencies are not parameters of this theory. 

3. Direct size characterization of multipolar SP modes 

In order to derive parameters of SP as a function of particle size, we use an accurate elec- 
trodynamic approach based on [49|], and described in more details in e.g. 50]. The formalism 
of Mie theory is used. However, the problem is formulated in absence of the illuminating light 
field; we are interested in intrinsic eigenproperties of the spherical particle (the analogy to the 
cavity eigenproblem) (3, 27]. In the present paper we extend the numerical calculations up to 
the radius of lOOOnm and plasmon polarity from I — 1 up to I — 10. That allowed us to explore 
new features of localized SP excitations and reconsider the results of . 

We consider continuity relations at the spherical metal/dielectric interface for the electro- 
magnetic (EM) fields, which are the solutions of the Helmholtz equation in spherical coordinates 
(r, 9, 0) inside and outside the sphere. The transverse magnetic (TM) modes of EM waves lo- 
calized on the interface possess a nonzero component of the electric field normal to the surface 
E r , which can couple with free-electron charge oscillations at the boundary r = R. The condi- 
tions for the nontrivial solutions of the continuity relations for TM mode define the dispersion 
relation: 



\Je in (u)^ (k out (u)R) ipi (k in (uj)R) - \js ut{u)£i {k out (uj)R) ip' t (k in (u)R) = (1) 

which is fulfilled for the complex eigenfrequencies of the field Qi in successive multipolar 
modes I = 1,2,3,... . ipi (z) and (z) are Riccati-Bessel spherical functions, the prime symbol 
(/) indicates differentiation with respect to the argument, k\ n = ^^Je in (uj) and k out = ^^s out , 
Einiuj) and e out are the dielectric functions of a metal and of the particle dielectric surrounding, 
respectively. The dispersion relation (EquationflTJ) is solved numerically for complex values of 
Qi{R) = u[{R) + iu'{(R) (u'{(R) < 0) for I in the range of 1 -=- 10. We solved the problem 
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for successive R values from 1 to lOOOnm in lnm steps. The f solve function of the MatLab 
program, utilizing the Trust-region dogleg algorithm was used. 

3.1. Material properties of plasmonic particles 

The simplest analytic function often used to describe a wavelength dependence of optical 

54 



properties of metals like gold or silver 
Sommerfeld model: 



51 



52 



53 



55| results from the Drude-Lorentz- 



e D (u) = e G - uj 2 p /{u 2 + i^u) (2) 

with e a = 1 (e.g. 

HQ). 

The effective parameter Eq > 1 describes contribution of bound 
electrons, u p is the effective bulk plasma frequency which is associated with effective concentra- 
tion of free-electrons, 7 is the phenomenological damping constant of electron motion. For bulk 
metals 7 = •fbuik is related to the electrical resistivity of the metal and is supposed to include all 
microscopic damping processes due to photons, phonons, impurities and electron-electron inter- 
actions. In the present modeling we accept the following effective parameters of the dielectric 
function of the bulk gold: e = 9, 84, u p = 9, OlOeV^, 7^ = 0, 072eV. 

The collision time l/^buik determines the electron mean free path in bulk metals. At room 
temperatures the electron mean free path in gold is 42nm [l|. It can be comparable or larger 
than a dimension of a particle. An additional relaxation term added Jo the relaxation rate 7 
accounts for the effect of scattering of free electrons by the surface [l 
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1r{R) = Ibuik + A 



R 



(3) 



where vf is the Fermi velocity (vf = 1.4 • 10 _6 m/s in gold), and A is the theory dependent 
quantity of the order of 1 [l|. We accept the value A = 0.33 in our modeling, according to 57]. 
Then the size-adopted bulk dielectric function is: 



e 7 ( fl )(w, R)=e - uI/{uj 2 + i-f R (R)u) (4) 

The real part of s 1 (r) is not modified by surface scattering. Frequency (and radius) depen- 
dence of the dielectric function (Equation [2] or HJ) couples to the overall frequency (and radius) 
dependence of the plasmon dispersion relation (Equation([T])). 

To underline the crucial role of interface damping in smaller particles, we will compare the 
results obtained with the dielectric function e in = ed{oj) given by Equation (J2J) (with 7 = jbuik, 
surface scattering neglected) and e in = e 7 ^m(u,R) given by Equation (j4]) (with 7 = j R (R)). 
The dielectric function of the particle surrounding (vacuum/air) is assumed to be e out = 1, or 
is chosen to reflect the index of refraction of the particle environment, as described in Section 



where we compare the results of our modeling with some experimental results of 34j, |37] 
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3.2. Resonance frequencies and damping rates vs radius; the role of the interface damping 

The solutions of the dispersion relation (EquationflTJ) define the size dependence of multipo- 
lar plasmon oscillation frequencies oj[(R) = Re(Qi(R)) and damping rates |u;"(.R)| = | Im(fli(R))\ 
(and SP damping times Ti(R) = h/\u"(R)\) of the EM surface modes oscillations. We express 
lj'i, uj'i, u and 7 in electronvolts (eV) for convenience. 

The size dependence of ojftR) and of SP modes for the first ten multipolar plasmons 

for gold spheres in vacuum/air (e out = 1, 7 = "/bulk = 0.072e\^, the radius up to lOOOnm) is 
presented in Figure [2 Black line (1 = 1) represents the dipole resonance frequency. Resonant 
excitation of SP oscillations takes place when the frequency of incoming light field u approaches 
the eigenfrequency of a plasmonic nanoantenna of a given radius R: u = u' t (R), I = 1,2, 3, .... If 
excited, plasmon oscillations are damped at corresponding rates (Figure [2b)). Damping 

of surface plasmon oscillations is the inherent property of SP modes that defines absorbing and 
scattering properties of the plasmonic particle as a function of size. 




R [nm] R [nm] 

Figure 2: a) Multipolar plasmon resonance frequencies <jj[(R) and b) plasmon damping rates |w"(i?)| (left axis) 
and corresponding damping times T\ (right axis) vs particle radius R. Reduction of the plasmon damping rates 
|w"(i?)| below the nonradiative limit 7/2) is demonstrated in the inset. 

The SP frequencies u[(R) decrease with size monotonically, as illustrated in Figure Ek). For 
given particle radius R, oj[(R) increase with plasmon polarity /. The multipolar SP resonances 
for larger particles are better spectrally resolved than those for the smaller ones. 

In the limit of small size, the analytic expressions for uj[(R) and ou"(R) can be found from 
the relation (JTJ) after applying the power series expansion of the spherical Bessel and Hankel 
functions. Keeping only the first terms of the power series one can get: 
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(6) 



For a perfect free-electron metal (eq — 1, 7 = 0) and e OU ( = 1, Equation dp leads to 
the well-known plasmon frequencies within the "quasistatic approximation" [l|, 0, 59|: ui' ol = 
u p / ^1 + Eout(l + 1)//, and in particular to the giant Mie resonance frequency u;q Z=1 = u p / y/3 
for I = 1. 

In the smallest particles, the plasmon damping rates |a; z "(i2)| fall to pure nonradiative damp- 



ing rates |wq, 



7 



nrad 



7/2, as demonstrated in Figure [2]d). The values of | c^q z I are the same 



for all multipolar plasmon modes I = 1,2,3... (Equation (jHJ)). The nonradiative damping rates 
are due to the ohmic losses if interface damping is neglected (7 = jbulk), or are equal to 



7 



nrad 



-fn/2 (Equation |3]) if interface damping is included in the modeling. 
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Figure 3: a) Dipole, quadrupole and hexapole plasmon resonance frequencies w[{R) and b) corresponding plas- 
mon damping rates |a>"(i?)| calculated without (solid lines) and with (lines with hollow circles) interface damping 
taken into account. Insets are magnifications of figures a) and b) for particles of small radii. 

Interface damping practically does not influence the size dependence of w[(R) as illustrated 
in Figure [3^,) . The minute red shift for the smallest sizes due to interface damping is s 
in inset of Figure [3k.). A similar finding results from the experimental observations 



nown 
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Therefore, the red shift of the (dipole) SP resonance frequency sometimes observed in small 
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particles versus decreasing size must result from some other phenomena due to the complex 
chemical effects and uncertainties in experimental samples 60 |. 

In extended size range that we study, plasmon damping rates are not simple mono- 

tonic functions of the radius (Figure Gb))- The initial fast increase of \u"(R) | with size is followed 
by gradual decrease for sufficiently large radii R. The increase in \ui +1 (R)\ in the subsequent 
I + 1 SP modes is followed by the decrease in the SP damping rates |w"(i?)| of lower polarity 
modes. 

The surface scattering effect affects the total damping rate |w"(-R)|, as illustrated in Figure 
[3b) and in the magnification of the part of this graph presented in Figure HI The surface 
scattering contribution to the total relaxation rate Jr(R) looses its importance for particles of 
radius larger than ~ 300nm (using the criterion (7^ — j)/jr ~ 1%). 
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Figure 4: Demonstration of the reduction effect of plasmon damping rates calculated without (7 = jbuik, 

solid lines) and with (7 = jr(R), lines with hollow circles) interface damping; the magnification of a part of 
Figure [31 



4. Radiative and nonradiative contributions to the total SP damping rates 



The processes leading to damping of SP have been the subject of wide discussion and ex- 
tensive studies, e.g. [l|, E2, 34, [fill, 62, 63] which, however, were limited to the dipole case. It 
was accepted that (dipole) SP in metal nanoparticles decays through both inelastic processes 
and elastic dephasing process which were usually neglected. Inelastic processes can be further 



divided into radiative and non-radiative decay processes [33j, |34|, [6J|. The SP damping time 
T was determined from the homogeneous linewidth r = 2h/T of a maximum in the spectrum 
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due to the SP dipole resonance. The observed spectral broadening of the maxima with size due 



to radiation is sometimes expected to be proportional to the particle volume [34j . 1571 . l64j |. Our 
results show (Figure [2b)), that the multipolar plasmon damping rates |u;"(i?)| are not simple, 
monotonic functions of particle radius. The initial fast increase of the total SP damping rate 
\u}'{(R)\ with size is followed by gradual decrease for sufficiently large radii for the given mode 
I. Our modeling suggest that the proportionality of the multipolar radiation damping rate to 
R 3 is not a general rule (see Figure [2b)). 

Our extended modeling allows us to consider the higher order multipolar damping rates 
|a;"(-R)| 1=1, 2. ..10, as well. The total multipolar damping times Ti/Ti can be related to the 



corresponding damping rates (homogeneous linewidths of the absorption [27J spectra Ti(R)) 
and their size dependencies in a natural way: 

V l {R) = 2\u'({R)\ = 2%/T l (7) 

and 

\ u i I = li +7 (8) 

where ^ nrad = r y/2. The assumption, that the nonradiative and radiative processes (in 
dip ole p lasmon mode) are additive and are independent, is commonly accepted (see for example: 



27 
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571 . l58l . l6lj|). However, expectation that the size dependence of the total damping rate 
\ui(R)\ is due to the size dependence of the radiative damping rate 7[ ad (i?) only (\u'{(R)\ = 
7[ ad (i?) + rynrad^ j^-j must be reconsidered (see Section [5] bellow). 

If the surface scattering effect is included, |o;f(i2)| becomes a decreasing function of size, 
starting from smaller particles (see Figure [3b) and HI lines with symbols). The total SP damping 
rates follow the size dependence of nonradiative damping: jr(R)/2: \u)"(R)\ ~ jr(R)/2 = 
•y nrad (R) in the range of radii which extends to larger R for growing SP polarity I. After 
reaching the minimum, |a; ; "(i?)| tends to follow the radius dependence unaffected by the interface 
damping. 

A contribution of the radiation damping ^[^(R) to the total damping rate |w"(i?)| can be 
treated as a measure of the ability of the particle to couple to the incoming field and to emit 
light in plasmonic mechanism. As long as |o;"(-R)| ~ ^nrad^ SP mode / can only weakly 
couple with the incoming radiation and has weak radiative abilities. Consequently, the weakly 
radiative plasmons appear in the absorption spectra with a smaller amplitude and in scattering 
spectra (see Figure [1]) they hardly manifest. With increasing I, SP plasmons gain the radiative 
character starting from larger particles due to the fact that 7™f(-R) > 7™f(i2). Therefore, the 
higher order SP modes can be excited (and appear in the scattering spectra) for larger particles. 
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This fact is known from Mie scattering theory, but its physical sense have not been explained 
by Mie solutions. 

5. Effect of reduction of multipolar plasmon damping rates 

In the easier case (interface damping omitted in the modeling), reduction of the total SP 
dumping rates j nrad bellow the 7/2 value (see the inset in Figure [2b)) manifests for plasmon 
modes with I > 1. This reduction effect can be clearly distinguished from other size dependent 
processes in some size ranges which grow with increasing I; for the quadrupole (/ = 2) plasmon 
mode, the reduction effect extends up to R ~ AOnm, for the hexapole (I = 3) plasmon mode up 
to R ~ 92nm, for / = 4 up to R ~ 155nm and so on. 

We have checked, the reduction of the total damping rate is not present if 7 = ([2]). On 
the contrary, it does not disappear if 7 7^ and Eq — 1. That suggests, that if the additive form 
of the Expression (JSj) holds, reduction of the total plasmon damping rate is due to the decrease 
of nonradiative decay j nrad as compared with its low size limit j nrad = <-y/2 resulting from 
the energy dissipation of meatal (absorption). Consequently, the radiative and nonradiative 
contributions have to be coupled by their size dependence. The model with interface damping 
included, reproduces the SP rate reduction below ■y nrad = •y R /2 (see Figure Hj), as well. 

We can conclude, that reduction of |o;f(-R)| below the 7/2 = jbuik/2 (or 7/2 = 7r/2) value 
takes place in the regions of sizes where the nonradiative damping is still not dominated by the 
fast radiative damping. Reduction in the total plasmon damping rate must be connected with 
suppression of the nonradiative decay channel. It can be inferred then, that the radiative and 
nonradiative processes are not independent. The effective rate of the nonradiative damping must 
be size dependent: ^ nrad = ^ nrad [R) with the value for the small particle limit: >y™ rad = 
Suppression of the nonradiative damping is not restricted to the size ranges for which the effect 
was demonstrated, but influences the optical properties of plasmonic particles in a large range 
of sizes. Reduction of 7™ raci (i?) manifests in the absorption spectrum of particles; absorptive 
abilities of large particles are poor as compared with small particles. This fact is described by 
solutions of Mie scattering theory, but its physical meaning have not been explained. 

The effect of the total damping rate reduction with particle dimension was observed for 
the first time in the experiment with gold nanorods [34|. The size dependence of the dipole 
plasmon damping rate was deduced from the homogeneous linewidth r = 2h/T of the maxima 
in the scattered intensities both for nanorods with various aspect ratios and spherical nanopar- 
ticles with radii in the range from 10 to 75nm. For nanorods, (dipole) plasmon damping rates 
decreased significantly for lower dipole plasmon oscillation frequencies. However, the decrease 
in the damping rate was not ascribed to the radiative processes. The authors of explain this 
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effect as reduction nonradiative plasmon decay by the fact that interband excitations in gold 
require a threshold energy of about 1.8eV and expect, that suppression of the damping rate 
in such mechanism would also be present in gold spheres, but for plasmon resonance energies 
below 1.8eV. Suppression of the plasmon damping rates in spheres had not been observed 
experimentally, as far as we know; the conclusive experimental data are limited to the plas- 
mon dipole only, while according to our study the effect manifests for / > 1. Importance of 
the threshold energy at 1, 8eV in gold and exclusion of the radiative damping as a reducing 
mechanism was not confirmed by our study. On the contrary, the reduction of multipolar SP 
nonradiative damping rates results from competition between radiative damping / y'[ ad (R) and 
all other damping processes included in 7 nrad (i?) and is present in plasmonic particles of any 
material. 



6. Quality factor of multipolar plasmon resonances 



Enhancement of optical response in resonance is usually described by the quality factor 
Q, defined as the product of the resonance center frequency and the bandwidth. In case of 
)lasmonic particles the quality factor is interpreted as a measure of the local field enhancement 
34j and is expected to define the effective susceptibility in nonlinear optical processes 36]. 
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Figure 5: a) Damping rates \oj['\ vs uj[ for successive radii R and b) quality factors Qi = uj[(R)/(2 |o;"(i?)|) for 
multipolar plasmon modes of gold nanoparticles (n ut = 1, surface scattering neglected). Graphs on the right 
are magnifications of circled regions of the graphs on the left. 



The dependence of the u'{ vs uj[ for successive radii R and the resulting quality factors 
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Qi(R) = oj[(R)/(2 \uj"(R)\) for multipolar plasmon modes of gold nanoparticles (n out = 1) are 
presented in Figures [5] (surface scattering neglected) and [6] (surface scattering included). The 
graphs on the right are the magnifications of the graphes (circled regions) on the left. 
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Figure 6: a) Damping rates \w'/\ vs u)[ for successive radii R and b) quality factors Qi = u)[{R)/{2 for 
multipolar plasmon modes of gold nanoparticles (n ut = 1, surface scattering included). Graphs on the right are 
magnifications of circled regions of the graphs on the left. 

Figure |5k) illustrates the effect of the total plasmon damping rate reduction to the value 
bellow ^nrad _ Q32e\Z for modes with I > 1. Figure Eb) illustrates the resulting quality factors 
Qi for successive SP modes which we use as a measure of energy stored in the SP oscillation of 
plasmon mode / at the resonance frequency oj[(R). In some ranges of SP resonance frequencies, 
the higher polarity plasmon modes (/ > 1) are more efficient in storing the SP oscillation energy. 
The similar conclusion comes from Figure Eb) (interface damping included). Quality factor of 
the dipole plasmon resonance reaches maximum value Qi=i ~ 29 for u>i{R) ~ 2, 606eV^ for gold 
nanoparticle of the radius R ~ lQnm (see Figure [7]) and is the fast decreasing function of size 
for both smaller and larger particles. If a gold particle is embedded in a medium of higher 
optical density (see Figure HH n out = 1,5), the maximum quality factor of the dipole plasmon 
resonance is even smaller, Qi=\ ~ 26. Such value corresponds to the resonance frequency of a 
gold particle of radius R ~ llnm. 
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Figure 7: Quality factors Qi vs radius R for multipolar plasmon modes of gold nanoparticles (n out — 1, surface 
scattering included). 

7. Dipole plasmon damping rates vs size: comparison with experimental results for 
silver and gold nanoparticles [34, 1371 ] 

Experimental investigations of spectral properties of plasmonic particles are usually per- 
formed on large particle ensembles, where inevitable variations in size, shape and surface prop- 
erties tend to mask the spectral properties of the individual particles. The effects of surface 
chemistry (especially in Ag particles) and uncertainties in sample parameters affect linewidths 



and maxima positions w 
chosen two experiments 



lich define damping times and resonance frequencies of SP. We have 



34 



37J, which provide the spectroscopic data for individual spherical 
particles of silver and gold. The experiments were performed in well-controlled particle environ- 
ments as a function of size for nanoparticles up to R = 25nm for silica-coated silver (Ag@SiC>2) 
and for gold nanoparticles immersed in a index matching fluid (n out = 1.5) for relatively large 
range of radii (up to R = 75nm). However, well resolved data was reported only for the dip ole 



plasmon resonance. In Figures |8] and M we compare the experimental results presented in 34, 3?J 
with the results of our modeling. 

Figure El illustrates the spectral width T measured in different single Ag@SiC>2 nanoparticles 
vs the inverse of their equivalent diameter D eq , optically determined by fitting the extinction 



spectra (see Figure 4 in |37j). Dashed line represents jr((2R) x ) dependence (Equation (j3J) for 



A = 0.7 and •jbuik = 0.125eV). These parameters [37J result from a linear fit to the experimental 



data. Solid line results from our modeling in the extended range of sizes for u p = 9.10eV, 
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Figure 8: Spectral width T of plasmon resonances measured [37] in single spherical Ag@SiC>2 particles vs the 
inverse of their equivalent diameter. Dashed line represents 7^((2i?) _1 ) dependence (equation ((S^-with the 
parameters A = 0.7 and "/bulk — 0.12561^ resulting from linear fit the line to the experimental data [371 ] . Solid 
line results from our model with the parameters uj p — 9.10, eo = ^eV, jbuik = 0.125eV and A = 0.7. 



e = 4eV, ■jbuik = 0.1256^ and A = 0.7. Our model with such input parameters reproduces the 
experimental data for smallest particles perfectly and describes the departure of cu/ = i((2i?) _1 ) 
dependence from linear. The oji(R) model functions describe consistently the size dependence 
of the SP damping rates for experimentally available particles and predicts damping rates (and 
the resulting damping times) for larger sizes (experimentally unavailable so for) and for higher 
plasmon multipolarities. 

Figure EK) illustrates spectral widths T vs resonance energy u' l=l (R) measured for different 
single Au nanoparticles (up to 2R = 150nm) (Figure 4 of 34 1). The experimental data seems 
to suggest that T decays linearly with the SP resonance frequency. However, it is not the case. 
The dependence of Ti = i(uii = i(R)) = 2\u" =1 (R)\ (line in figure Ek)) reproduces the experimental 
data and predicts the decrease in Ti =1 for particles larger than those studied experimentally. 
Such decrease in damping rate |o; / " =1 (i?)| of the dipole SP is accompanied by the increase in 
the damping rate of the quadrupole \cui =2 (R)\ SP and of the following higher polarity SPs, as 
discussed in Section [5] (see Figure Eb))- 

Figure [9b) illustrates the quality factor of the dipole resonance Qi=\ = u' l=1 (R)/2\u" =l (R)\ 
vs resonance energy u' l=l (R) (solid lines). Our modeling shows that interface damping (dashed 
line) substantially reduces the quality factor for particles from the smallest size range. The 
maximum factor Qi=± ~ 26 is found in nanoparticles of radii R = 12nm. The agreement with 
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the experimental data of 



34 ] is very good. 




Figure 9: a) Spectral widths T of plasmon resonances vs resonance energy measured [34j | in single spheri- 
cal Au particles, b) Quality factor resulting from the experimental data. Lines represent the spectral width 
r/ = i(wi = i(i?)) = 2\uj" =1 (R)\ and the quality factor Qi—% = u)' l _ l (R)/2\ui'/_ 1 (R)\ vs lj' 1=1 (R) for the dipole mode 
obtained from our modeling with surface damping neglected (dashed line) and included (solid line). 



8. Conclusions 

The dependence of the SP resonance frequencies wftR) and damping rates \cui(R) \ on particle 
radius determine the intrinsic optical properties of plasmonic spheres (illuminated or not). Our 
study provides direct, accurate size characteristics for a broad range of particle radii (up to 
R = lOOOnm) and plasmon polarities (up to I = 10). At present, this exceeds the range 
experimentally explored. 

Size dependence of SP damping rates |u^"(i2)| allows to distinguish the size ranges in which 
efficient transfer of radiation energy into heat takes place (large contribution of the nonradiative 
decay) and those in which the particles are effective radiating antennas (dominant contribution 
of the radiative damping). As long as the contribution of the radiative decay is negligible, the 
particle is not able to couple to the incoming field effectively and has weak radiative abilities. 
If the contribution of the radiative damping prevails, the particle is able to emit light within 
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plasmonic mechanism efficiently. Effective radiating enables efficient interaction of SP near- 
field with other structures at a desired resonance frequency u = ojI(R). The knowledge of 
size dependence of both: the multipolar SP resonance frequencies u)[(R) and the corresponding 
damping rates \ui(R)\ = h/Ti(R) is indispensable to shape the particle plasmonic features 
effectively 

Our study, extended toward large particle sizes and plasmon multipolarities, revealed new 
features of the total plasmon damping rates. In certain ranges of radii, reduction of the multi- 
polar SP damping rates, as compared with its small size limit, takes place. The small size limit 
is equal to the nonradiative damping rate resulting from absorption and heat dissipation. The 
suppression of nonradiative damping is not limited to the size range, in which the effect was 
directly demonstrated. It is present when the radiative damping brings the dominant contribu- 
tion to the total plasmon damping. The reduction of ■y nrad (R) with the growing contribution of 
radiation damping is revealed in the absorption spectra of particles; absorptive abilities of large 
particles are poor. This fact is described by solutions of Mie scattering theory, but its physical 
background have not been explained, as far as we know. Our study led us to conclusion, that 
the reduction of multipolar SP nonradiative damping rates results from competition between 
radiative damping and all other damping processes included in 7 nrad (i?). As far as we know, 
such hypothesis has not been proposed before. Our study provides a new starting point for 
better understanding of rules that define contributions of plasmon radiative and nonradiative 
decay rates to the total SP decay rate as a function of particle size. The independent modeling is 
necessary for better understanding of the role of SP radiative and nonradiative decay channels. 

The proposed SP characteristics vs particle size provide a consistent, uniform description 
of the experimental SP damping rates measured for single spherical particles in [3^] and 



37] . 



Not only dipole damping rates (found experimentally) but also multipole damping rates (and 
resulting damping times) in a broad range of sizes can be predicted. 

The derived size dependence of plasmon decay rates |u;"(i2)| at resonance frequencies uj[(R 
), and the quality factors Qi(R) of SP multipolar modes define not only the spectral scattering 
abilities of the plasmonic spheres, but also reflect changes in the strength of coupling of SP 
modes with the external field. Such characteristics can serve as a tool for controlling the spectral 
features of plasmonic nanospheres in technological or diagnostic applications, by optimizing the 
size at the desired resonance frequency u = oj[(R). 
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